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Low-lying relaxation satellites give rise to ultrafast electron dynamics following outer-valence 
ionization of a molecular system is studied. To demonstrate the impact of such satellites, the 
evolution of the electronic cloud after sudden removal of an electron from the highest occupied 
molecular orbital (HOMO) of the organic unsaturated nitroso compound 2-Nitroso[l,3]oxazolo[5,4- 
d] [l,3]oxazole is traced in real time and space using ab initio methods only. Our results show that 
the initially created hole charge remains stationary but on top of it the system reacts by an ultrafast 
7r — 7r* excitation followed by a cyclic excitation-deexcitation process which leads to a redistribution 
of the charge. The n — n* excitation following the removal of the HOMO electron takes place on 
a sub-femtosecond time scale and the period of the excitation-deexcitation alternations is about 
1.4 fs. In real space the processes of excitation and de-excitation represent ultrafast delocalization 
and localization of the charge. The results are analyzed by a simple two- and three-state model. 



I. INTRODUCTION 



With the advent of the attosecond pulse techniques 
(see, e.g. Refs. [l|, [|| and the references therein) the sci- 
entific community obtained a powerful tool to monitor 
the electron dynamics in atomic and molecular systems 
and to study processes that take place on a time scale 
in which the electronic motion is still disentangled from 
the slower nuclear dynamics. Such kind of processes are, 
for example, the response of the electronic cloud to an 
ultrafastperturbation, like ultrafast excitation or ioniza- 
tion 

I, H,H 10- About ten years ago in a first work 
Q it was shown that indeed after a sudden ionization 
rich ultrafast electron dynamics may occur. The positive 
charge created after the ionization can migrate through- 
out the system solely driven by the electronic many-body 
effects - the electron correlation and electron relaxation 
[1, Q. Since this charge mi grat ion is ultrafast, typically 
few femtoseconds 0, E3, [TTI Il2l. Il3j , it can be calculated 
neglecting the nuclear motion as long as one is concerned 
with the relevant time interval during which this ultrafast 
process takes place. Clearly, at later times the nuclear 
dynamics will come into play and will perturb the pic- 
ture. Thus, if we wish to know precisely what happens at 
this later stage, the nuclear motion must be considered. 
However, since several or even many electronic states par- 
ticipate, an adequate description of the nuclear motion is 
rather involved, ft should be noted, in this respect, that 
very recently a time-dependent Born-Oppenheimer ap- 
proximation for treating quantum mechanically coupled 
electron-nuclear dynamics was proposed [141 . Ila ] which 
might provide a possibility to attack such complicated 
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problems. 

We have studied the charge migration in different 
molecules [|, ES El El El and found that it usu- 
ally takes place after ionization in the inner-valence shell 
where the electron correlation effects are typically much 
stronger than in the outer-valence shell. However, very 
recently we have shown that the charge migration phe- 
nomenon is not inherent only to the inner- valence ionized 
states but can take place also after an outer-valence ion- 
ization pjj]. The system studied in Ref. [l6| possesses a 
chromophore-donor site which is initially ionized and the 
created hole charge migrates throughout the system to 
the amine-acceptor site within just 4 fs. Further analysis 
[ijj then showed that the mechanism underlying this ul- 
trafast migration is the so-called hole mixing [l8| which 
is one of several ways of the manifestation of electronic 
many-body effects in the ionization process 

In the present work we want to continue the investiga- 
tion of ultrafast electron dynamics following ionization of 
the outer-valence shell focusing on a different mechanism, 
namely the so-called dominant satellite mechanism. It is 
known from the early days of the photoelectron spec- 
troscopy that the correlation and relaxation effects can 
lead to the appearance of additional weak bands in the 
photoelectron spectrum, the so-called shake-up or satel- 
lite bands, which correspond to excitation processes ac- 
companying the ionization 20]. This kind of shake-up 
states are typical for the inner-valence and core ioniza- 
tion and appear in the ionization spectra of nearly ev- 
ery atom or molecule. However, numerous theoretical 
and experimental studies showed that there are certain 
classes of compounds where satellites can appear also in 
the outer-valence region. These are systems containing 
heavier atoms, like transition-metal complexes j2ll.l22j]. or 
7r-electron systems @, H3, [H, HI H3 , which have low- 
energy virtual orbitals. Here we will concentrate on the 
latter class of systems and will investigate the impact of 
the low-lying satellites on the electron dynamics follow- 
ing outer- valence ionization of some organic unsaturated 
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nitroso compounds [27| . 

The paper is organized as follows. In Sec. |TT] the 
theoretical background of the methodology used for 
calculating the ultrafast electron dynamics is briefly 
outlined (Sec. Ill A[) and the different mechanisms of 
charge migration are briefly discussed (Sec. Ill Bl) . 
In Sec. Mil we present the results of our calcula- 
tions including the ionization spectrum of the molecule 
2-Nitroso[l,3]oxazolo[5,4-d][l,3]oxazole and the electron 
dynamics following ionization of its highest occupied 
molecular orbital (HOMO). Section ITVl is devoted to the 
discussion and analysis of the results obtained via a sim- 
ple two- and three-state model. 

II. THEORETICAL BACKGROUND 
A. Basic equations 

In this section we briefly review the theoretical back- 
ground of the methodology used to study ultrafast elec- 
tron dynamics following ionization of a system. For tech- 
nical details we refer the reader to Refs 0, [ll|, [H[ . 

The starting point of our investigation is a neutral 
molecule in its ground state \^o}- The ionization of 
the system brings it into a non-stationary state \<&i). A 
convenient quantity then for tracing the succeeding elec- 
tron dynamics is the density of the so created initial hole 
which can be defined by the following expression: 

Q(r,t) := {y \p(r,t)\* )-{$iW,t)\®i) = Poff-piftt), 

where p is the electron density operator. The first term 
on the right-hand side of Eq. JT]), po, is the ground state 
density of the neutral system and is time-independent. 
The second term, pi, is the density of the cation and 
hence is time dependent, since is not an eigenstate 
of the cationic system. The quantity Q(f,t), referred 
hereafter as the hole density, describes the density of the 
hole at position r and time t and by construction is nor- 
malized at all times t. 

In the Heisenberg picture, the time-dependent part 
Pi (r , t) reads: 

Pl (r,t) = {^\e i6t p(f,0)e- i6t \^) = (^(t)WMHt)), 

(2) 

where | $»(*)) = e~ iHt \^i) is the propagating multielec- 
tron wavepacket. 

Using the standard representation of the density oper- 
ator in a one particle basis tp p (r) and occupation numbers 
n p , Eq. {T]) can be rewritten as follows 

Q(r,t) = J2^ P (^ q (r)N Pq (t)- (3) 

where the matrix N(i) = {N pq (t)} with elements 

N pq (t) = 6 J2(Mt)\^M)pMN(^N\Mt)) (4) 

M,N 



is referred to as the hole density matrix. The second term 
of Eq. ([!]) is obtained by inserting in Eq. © a resolution 
of identity of a complete set of appropriate ionic eigen- 
states I^m) before and after the density operator p(r, 0). 
The matrix pMN is thus the representation of the density 
operator within this basis. 

Diagonalization of the matrix N(t) for fixed time 
points t yields the following expression for the hole den- 
sity 

Q(r,t) =^|^ p (r,i)| 2 n p (i), (5) 

p 

where <p p (r, t) are called natural charge orbitals, and h p (t) 
are their hole occupation numbers. The hole occupation 
number h p (t) contains the information on which part of 
the created hole charge is in the natural charge orbital 
tpp(f,t) at time t. Because of the charge conservation, 
one finds that ^ p n p (t) — 1 at any time t. 

For calculating the hole density and its constituents 
we use ah initio methods only. The whole calculation 
consists of four steps. After determining the molecular 
geometry, the first step is a Hartree-Fock (HF) calcula- 
tion. The second step is the calculation of the relevant 
part of the ionization spectrum via Green's function for- 
malism. A computationally very successful approach to 
obtain the Green's function is the algebraic diagrammatic 
construction [ADC(n)] scheme [28j . In the present cal- 
culation we used the non-Dyson ADC(3) method [29l.[30j 
realized within the so-called intermediate-state represen- 
tation [3l|, [HJ , an effective many-body basis serving as 
I^m) introduced in Eq. ([J|. The third step is the prop- 
agation of the multielectron wavepacket of the ionized 
system [ll[ with the help of the short iterative Lanc- 
zos technique (33|. The fourth and last step is to build 
through Eq. (01 the matrix N(f) and to diagonalize it in 
order to obtain the natural charge orbitals (p p (r , t) and 
the hole occupation numbers n p (t), see Eq. ([5]). With the 
help of these quantities we can now trace the evolution 
of the hole density of a system after suddenly removing 
one of its electrons. 

At this point we would like to comment briefly on 
the choice of the initial state The above sketched 

methodology is independent of the particular choice and 
the way of preparation of the initial state as long as the 
ionized electron is removed from the system on a shorter 
time scale than that of the charge migration. The as- 
sumption made is that the initially created ionic state 
can be described by a separable manyelectron wavefunc- 
tion, i.e. the interaction between the ionized electron 
and the remaining core is neglected - sudden approx- 
imation (see, e.g. [IH). Within this approximation 
the initial hole is described by the so-called Dyson or- 
bital, i.e. the overlap between the iV-electron initial and 
(N — l)-electron final wavefunctions. However, in the 
outer-valence region the Dyson orbitals differ very little 
from the canonical Hartree-Fock orbitals [H, [3(| . That 
is why, in the numerical calculations to be discussed in 
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this paper the initial state is prepared through a sud- 
den removal of an electron from a particular HF-orbital. 
The initial state can, of course, be constructed such that 
it corresponds to a removal of an electron from a lin- 
ear combination of HF-orbitals. This liberty allows one 
to reproduce the hole density of practically every par- 
ticular initial vacancy. However, to avoid investigating 
many linear combinations of HF-orbitals of interest we 
concentrate on specific HF-orbitals. In this way we can 
unambiguously identify the basic mechanisms leading to 
charge migration. Since the time-dependent Schrodinger 
equation which governs the electron dynamics is a lin- 
ear equation, these mechanisms are also operative when 
other choices of the initial state are used. 



B. Ionic states and mechanisms of charge migration 

Here we will briefly review the basic mechanisms of 
charge migration referring the interested reader to Ref. 
0| for more details. 

For proper understanding of the basic mechanisms of 
ultrafast charge migration following ionization it is illu- 
minative to analyze a typical ionization spectrum. The 
calculated cationic spectrum consists of vertical lines, 
where each line represents a cationic eigenstate \I). The 
position of the line is given by the ionization energy, 
and its hight - by the square of the transition amplitude 
a quantity related to the ionization cross section. 
For the ease of interpretation we will expand the exact 
cationic state \I) in a series of electronic configurations, 
as is traditionally done in configuration interaction (CI) 
calculations (see, e.g. Ref. H3|): 

|/> = E c i° a il*°>+ £ 4h^iI*o> + ---, (6) 

j a,k<l 

where \^q) is the exact ground state of the neutral sys- 
tem, a and a) are the annihilation and creation oper- 
ators, respectively, and c'^'s are the expansion coeffi- 
cients. The indices a, 6, . . . refer to unoccupied (virtual) 
orbitals (or particles), whereas the indices i,j, ... indi- 
cate occupied orbitals (or holes). Throughout the whole 
text p,q, . . . will be used as general indices. Accordingly, 
the terms fijl^o) are called one-hole (lh) configurations, 
since one electron has been removed from the correspond- 
ing occupied orbital, the terms djafca/ l^o) are referred to 
as two-hole-one-particle (2hlp) configurations, indicating 
that in addition to the removal of one electron another 
one is excited to a virtual orbital, and so forth. Note that 
in the spirit of the GF approach, which accounts also for 
the ground-state correlations, in Eq. ^ the expansion 
is applied on the exact ground state l^o)) rather on the 
uncorrelated HF one |$o) as in the usual CI calculations. 

As was noted above, in our study the initially pre- 
pared non-stationary ionic state is created by suddenly 
removing an electron from a particular HF-orbital, i.e. 
by acting with the corresponding annihilation operator 



on the ground state of the neutral, = fijl'&o)- Thus, 
only the lh configurations contribute to the transition 
amplitude i.e. to the spectral intensity. Without 

correlation effects, the spectrum will consist of lines, one 
for every occupied orbital (pi, with intensities equal to 1. 
If correlation effects are weak, the ionization spectrum 
will consist of main lines which have large overlap with 
the lh configurations. This is typical when ionizing the 
outer-valence shell of a system. In this case the molecular 
orbital picture is still valid. If the correlation effects are 
stronger, beside the main line satellite lines will appear. 
The intensities of the satellite lines are weaker than those 
of the main lines since they correspond to cationic states 
that are dominated by 2hlp configurations and have only 
small or moderate overlap with the lh configurations. 
Three types of satellites can be distinguished |19j: re- 
laxation satellites where at least one of the two holes 
in the 2hlp configuration is identical to the lh orbital 
of the main line, correlation satellites where both holes 
differ from the lh orbital of the main line, and ground- 
state- correlation satellites stemming from the correlation 
effects present in the ground state of the neutral. In the 
inner- valence, where the correlation effects are strong, 
the distinction between the main lines and the satellites 
cease to exist and the spectrum becomes a quasicontin- 
uum of lines with small to moderate intensities. This 
phenomenon is known as breakdown of the molecular or- 
bital picture [l9| . 

Depending on the structure of the ionic states involved, 
three basic mechanisms of charge migration have been 
identified 0: the hole mixing case, the dominant satellite 
case, and the breakdown of the molecular orbital picture 
case. In what follows we will briefly describe them. 

(i) Hole mixing case. For simplicity, we will consider 
the two-hole mixing, i.e. the situation when two lines in 
the spectrum correspond to ionic states which are linear 
combinations of two lh configurations o^l^o) and a,k\^o)- 
In this idealized case, if we create the initial hole in one of 
the orbitals, say (pj, then the hole will oscillate between 
the two orbitals ipj and ip k with a frequency determined 
by the energy difference between the two ionic states. If 
the two orbitals are localized on two different sites of the 
system, the hole mixing mechanism will lead to an oscilla- 
tion of the initially created positive charge between these 
two sites. The hole mixing mechanism was identified as 
the driving force of the ultrafast charge migration follow- 
ing outer- and inner-valence ionization in many different 
molecular systems (see, e.g. Refs. Elm, HO). 

(ii) Dominant satellite case. Let us consider the sit- 
uation when we have two ionic states, a main state and 
a satellite, both having overlap with the original lh con- 
figuration. In CI language these states can be written 
as 

\I m ) = ci a,i\^o) + C2 ala k ai\^o), 
\I S ) = c 2 Oi|*o) - ci ala k ai\^ ), 

where the two coefficients C\ and c 2 satisfy the equation 
c\ + c\ = 1. In this idealized case, assuming that all in- 



volved orbitals ifi, (fk, tpi, and fa are different, i.e. the 
case of a correlation satellite, we will observe the fol- 
lowing electron dynamics succeeding the ionization out 
of orbital (fi. The hole initially localized on orbital tpi 
will migrate to the orbital tpk ( or fl) accompanied by 
an excitation from orbital tpi (or (fk) to the virtual or- 
bital tp a - Again the dynamics will be oscillatory with a 
period determined by the energy differences between the 
states \I m ) and |J S ). This type of mechanism was identi- 
fied to be responsible for the ultrafast electron dynamics 
taking place after inner-valence ionization in several sys- 
tems [ill H2I . |38| . We note at this point that there is also 
a dominant relaxation satellite mechanism which is not 
discussed so far and will be investigated for the first time 
in the present paper. 

(Hi) Breakdown of the molecular orbital picture case. 
In the inner-valence region of the spectrum, where the 
quasicontinuum of lines appears, one can distinguish two 
general cases depending on whether the states are be- 
low or above the double ionization threshold of the sys- 
tem. Supposing that the quasicontinuum of states has 
a Lorentzian shape, in both cases the initially ionized 
orbital will "lose" its positive charge exponentially with 
time. If the states are below the double ionization thresh- 
old, the charge will be typically shared among many other 
orbitals and at the end of the process will be spread more- 
or-less uniformly over the whole cation. This situation 
was studied in Refs. d, HI]. If the states are above 
the double ionization threshold, i.e. an electronic decay 
channel is open, this mechanism will describe the process 
of emission of a secondary electron. This situation was 
studied in Ref. [39j . 

In the present work we will concentrate on the domi- 
nant relaxation satellite mechanism of charge migration, 
and will discuss for the first time the manifestation of this 
mechanism after ionization of the outer-valence shell. 



III. RESULTS 

We applied the methodology sketched in Sec. Ill Al to 
several molecules known to posses low-lying satellites in 
order to trace in real time and space the response of 
the electronic cloud of these systems to the sudden re- 
moval of an electron from their outer-valence shell. The 
systems studied were organic unsaturated nitroso com- 
pounds theoretically investigated by Wardermann and 
von Niessen [27| ■ The molecular geometries were opti- 
mized using DFT methodology (BP86/SV(P)). Through- 
out the remaining calculations standard DZ basis sets 
[40l ] were used. All studied molecules show very simi- 
lar behavior, thus in this paper we will present only one 
of them, namely the molecule 2-Nitroso[l,3]oxazolo[5,4- 
d][l,3]oxazole (see the sketch). 

The molecule is planar and belongs to the C s sym- 
metry point group having 29 occupied orbitals in the a 1 
irreducible representation (irrep) and 6 occupied orbitals 
in the irrep a" . The computed outer- valence ionization 
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spectrum of the molecule is shown in Fig. [TJ The calcu- 
lations were performed via Green's functions (GF) based 
approach, namely the non-Dyson ADC (3) method [29j |. 
The states belonging to the a" irrep are plotted in black. 
All other colors are related to different lh-configurations 
resulting from the removal of an electron out of a spe- 
cific HF-orbital belonging to the a' irrep. It is seen that 
the main line at 9.6 eV, the satellite at 11.2 eV and the 
larger contribution to the satellite at 12.6 eV come from 
the same lh configuration, namely the 29a' _1 configura- 
tion which is given in green in the figure. The orbital 29a' 
is the highest occupied orbital (HOMO) of the molecule. 

1 | , , , , 1 
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FIG. 1: Outer- valence part of the ionization spectrum of 
the molecule 2-Nitroso[l,3]oxazolo[5,4-d][l,3]oxazole calcu- 
lated via the non-Dyson ADC(3) Green's function method. 
Each vertical line shown is related to a final cationic state 
and is located at the corresponding ionization energy. The 
states belonging to the a" irrep are plotted in black. All 
other colors are related to different lh-configurations result- 
ing from the removal of an electron out of a specific HF-orbital 
belonging to the a' irrep. 



Let us take a closer look at the cationic states corre- 
sponding to these three lines in the spectrum, the main 
line and the two satellites. The first cationic state, cor- 
responding to the line at 9.6 eV, is constructed mainly 
from the 29a'~ 1 lh configuration (63%) and from the 
29a'" 1 %a"~ 1 7a" 2hlp configuration (25%). The second 
state, corresponding to the satellite line at 11.2 eV, is also 
mainly constructed from the 29a' -1 and 29a' _1 6a" _1 7a" 
configurations but their weights are inverted compared 
to the first cationic state. The third state at 12.6 eV 
consists of 17% 29a'" 1 , 10% of 27a'" 1 , and 43% of 
29a'~ 1 6a"~ 1 7a". In all states the missing contributions 
are distributed over many other lh and 2hlp configura- 
tions. 
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Obviously, we encounter the situation of a main line 
stemming from the ionization out of the HOMO and two 
relaxation satellites corresponding to an excitation of an 
electron from HOMO-1 to the lowest unoccupied molec- 
ular orbital (LUMO). The three involved HF-orbitals 
(the HOMO, orbital 29a', the HOMO-1, orbital 6a", the 
LUMO, orbital 7a") are displayed in Fig. 03 It is seen 
that the HOMO is localized on the nitroso-site (N=0 
site) of the molecule, while the HOMO-1 and LUMO are 
delocalized ir orbitals. 

a) b) c) 




FIG. 2: Hartree-Fock orbitals of the molecule 2- 
Nitroso[l,3]oxazolo[5,4-d] [l,3]oxazole: a) orbital 29a', or the 
HOMO; b) orbital 6a", or the HOMO-1; and c) orbital 7a", 
or the LUMO. 

Let us now see what happens after a sudden removal 
of an electron out of the HOMO. In order to analyze the 
electron dynamics following the ionization it is illumi- 
nating to trace the time evolution of the hole occupation 
numbers n p (t). In Fig. [3] we present the evolution of the 
hole occupation numbers during the first 5 fs after the 
ionization. The occupation numbers corresponding to or- 
bitals belonging to irrep a' are shown in black, while those 
corresponding to orbitals belonging to a" in red. Only 
five of the occupation numbers from each irrep (those 
who contribute the most) are shown. At the beginning 
of the process all occupation numbers are equal to zero 
except one which corresponds to the initially ionized nat- 
ural charge orbital (pi, indicated by V in the figure. At 
that time the natural charge orbital <pt has 100% overlap 
with the HF-orbital 29a'. We see that apart from the 
very fast drop during the first 50 attoseconds, which is 
a universal response of a system upon sudden ionization 
(see Ref. 41]), the initially created hole charge stays 
in orbital ip>i throughout the whole studied period. One 
should keep in mind that the natural charge orbitals also 
change with time, but in this case the orbital (pi varies 
very little and at any time point has more than 95% 
overlap with the HF-HOMO. Thus, the positive charge 
practically stays where it is initially created. 

However, another process takes place. We see that for 
approximately 0.7 fs more than 80% of another hole ap- 
pears on the natural charge orbital denoted by "fc" in 
Fig. [3j and in the same time more than 80% of an elec- 
tron is promoted to the natural charge orbital denoted by 
"a" in the figure (note that "negative" hole occupation 
number describes a particle, i.e. an electron in a vir- 
tual orbital). The analysis shows that at all time points 
the natural charge orbital <pk overlaps more than 80% 
with HF-orbital 6a", while the natural charge orbital <pt 
overlaps more than 60% with HF-orbital 7a". Thus, the 
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FIG. 3: Hole occupation numbers as a function of time for 
the first 5 fs after ionization of the HOMO (orbital 29a') of 
the molecule 2-Nitroso[l,3]oxazolo[5,4-d] [l,3]oxazole. The five 
most important occupation numbers of each irreducible rep- 
resentation are shown. At t = all occupation numbers are 
equal to except of the one which relates to the initially 
ionized orbital. The negative hole occupations have to be in- 
terpreted as the corresponding fraction of an electron filling 
the respective virtual orbital. 



removal of an electron from the HF-HOMO leads to an 
ultrafast (less than a femtosecond) ir — ► it* excitation 
from orbital 6a" to orbital 7a". From Fig. [3] we see also 
that this excitation process is cyclic, i.e. one observes al- 
ternating excitations and de-excitations with two beating 
periods one nearly twice as long as the other. 

Before analyzing in more detail the mechanism under- 
lying this ultrafast electron dynamics, let us see how the 
process develops in space. For that purpose in Fig. [f] 
we display snapshots of the evolution of the hole den- 
sity, Eq. ([5]), at different time points covering the time 
of a full oscillation cycle. The points chosen are t = 0, 
0.7, 1.4, 2.1, and 2.8 fs which correspond to the minima 
and maxima of the curves k and a in Fig. [3J We see 
that at t = the charge is localized on the nitroso side 
of the molecule (N=0 group). As time proceeds a sec- 
ond hole starts to open at the left oxazole ring and an 
electron starts to appear at the nitroso group (shown in 
orange in Fig. [f]) . One has to keep in mind that since the 
orbitals 6a" and 7a" are delocalized at the places of over- 
lap we have a mutual cancellation between the electronic 
and hole densities. Thus, the process represents alternat- 
ing ultrafast derealization and localization of the charge. 
However, due to the large number of electronic states in- 
volved these alternations are not purely repetitive. From 
Fig. [|]it is clear that the charge tends to delocalize more 
and more as time proceeds. 
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FIG. 4: 3D hole density Q(f,t) at times t = 0, 0.7, 
1.4, 2.1, and 2.8 fs after ionization of the HOMO of 2- 
Nitroso[l,3]oxazolo[5,4-d] [l,3]oxazole. The "negative" hole 
density, or the electron density, is shown in orange. 

IV. ANALYSIS AND DISCUSSION 

Here we would like to analyze in more detail the 
underlying mechanism of the ultrafast electron dynam- 
ics observed after sudden ionization of the HOMO of 
2-Nitroso[l,3]oxazolo[5,4-d][l,3]oxazole. The dominant 
satellite mechanism studied in Ref. Q and described 
briefly in Sec. Ill B I provides an understanding for the im- 
pact of a correlation satellite on the electron dynamics 
following ionization. It leads to the migration of the ini- 
tially created hole charge and a simultaneous excitation 
on top of it. However, as we saw in the previous section, 
in the case of dominant relaxation satellite we observe 
only the excitation while the initially created charge re- 
mains stationary (see the red and black curves, denoted 
as k, a, and i, respectively in Fig. [3J). This behavior can 
be understand with the help of a simple model. 

Suppose we have only two cationic states - a main state 
and a relaxation satellite. In CI language these states can 
be written as 

\I m ) =Cl|*i) +C 2 \^aik), 

\I S ) = C 2 |*i) -ci|* oife ), 

where = a;|* ) and |* Qlfe ) = ata 4 a fe |<F ). After 
some trivial algebra (see Ref. one arrives at the fol- 
lowing analytical expression for the hole density matrix 
[Eq. ©] 

N(t) = 

/ 1 \ 

-2(cic 2 ) 2 [l - cos(wt)] -cic 2 (ci - c§)[l - cos(wt)] , 
\0 -eic 2 (cj - c|)[l - cos(wt)] 2(cic 2 ) 2 [l - cos(u)t)] / 

(7) 

where uj = (Ei s — Ei m ) / h. 

The diagonalization of this matrix gives the time- 
dependent hole occupation numbers (using that c\-\-c\ = 
1): 

hi(t) = 1, 

n k/a (t) =±2cic 2 sin 2 (^/2). (8) 



We see that indeed the initially created hole is station- 
ary and the occupations of the two other natural charge 
orbitals involved in the 2hlp configuration oscillate such 
that the total charge is equal to 1 at any time. 
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FIG. 5: Model time-dependent occupation numbers of the 
three natural charge orbitals (i, k, a) in the case of a domi- 
nant relaxation satellite. The initially created hole (i) remains 
stationary, while an excitation from the occupied orbital (k) 
to the virtual orbital (a) is observed. The following param- 
eters were used in the model: ci = V0.7, = \/0.3, and 
lo = 2.4 fs" 1 [see Eq. ©]. 



The time evolution of the three hole occupation num- 
bers in this model are shown in Fig. In order to com- 
pare with the full calculation presented in Fig. [3J the 
parameters used in the model are c\ = V0.7, c 2 = V0.3, 
and uj = 2.4 fs -1 , which corresponds to the situation re- 
alized in the studied molecule taking only the main state 
at 9.6 eV and the first relaxation satellite at 11.2 eV. It 
is seen that this simple two-state model reproduces the 
overall behavior of the cyclic excitation and de-excitation 
with a period of about 2.6 fs. The more involved beating 
pattern observed in Fig. [3] can be explained easily by the 
influence of the third state at 12.6 eV. 

It is straightforward to show that in the case of a main 
state and two relaxation satellites having the form 

iu=4 m) i* i >+4 m) i^>+4 m) i^>, 
i/ sl ) = 4 sl) i* i )+4 sl) i^)+4 sl) i* aife ), 

|/ s2 ) = c[ s2) \^) + 4 s2) i*,) +4 s2 

and supposing that S™^ = 4 — 0, i.e. the lh config- 
uration does not contribute to the main state and 
one of the two satellites, the hole occupation numbers 
describing the excitation k — > a will be proportional to 
[sin 2 (tJii/2) + sin 2 (cj 2 t/2)], where uj 2 = (E Isl - E Im )/h 
and u)\ = {Ei b2 — Ei m )/h, while the orbital j will get no 
occupancy, i.e. hj(t) = 0. This is exactly the case real- 
ized in the example studied in the present paper where 
the lh configuration corresponds to 27a' -1 which 
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contributes only to the second satellite state at 12.6 eV. 
It is seen that this second relaxation satellite will then 
just introduce an additional frequency 102 = 4.6 fs -1 cor- 
responding to an oscillation period of about 1.4 fs and 
we will get the beating oscillation behavior observed in 

Fig. m 

The following comment is in order. Although the relax- 
ation satellite does not affect the dynamics of the initial 
hole charge itself (the initial hole remains stationary), 
depending on the spatial distribution of the involved or- 
bitals the dominant relaxation satellite mechanism can 
lead to a charge migration via hole screening. If the elec- 
tron is ejected from a localized orbital i and the unoccu- 
pied orbital a is localized in the same region of space, then 
the initially created hole will be screened by the excita- 
tion k — > a. Thus, in the idealized case when the orbitals 
i and a are localized on one site of the system, while the 
orbital k is localized on a different site, the dominant re- 
laxation satellite mechanism will lead to oscillations of 
the hole charge between these two moieties of the sys- 
tem. It has to be noted that the unoccupied orbitals 
are often delocalized and the realization of this simple 
mechanism of charge migration from one site to another 
is not very likely in such systems. This is also seen in the 
example presented in this work where the process repre- 
sents rather an ultrafast derealization and localization 
of the charge. However, in larger and/or strongly corre- 
lated systems where many orbitals can be mixed in the 
ionic states the hole-screening charge migration mecha- 
nism may be operative. We note that there do exist many 
systems with localized unoccupied orbitals, e.g. para- 
Nitroaniline 0, H^|. Furthermore, in large systems even 
a* orbitals are likely to be located on individual moieties. 
In all these kinds of systems one could expect ultrafast 
charge migration from one part of the system to another 
via the hole-screening excitation mechanism. 



Let us summarize. In the present work we studied 
the ultrafast electron dynamics following outer-valence 
ionization of a molecular system concentrating on the 
impact of low-lying relaxation satellites. For that pur- 
pose we traced the evolution of the electronic cloud af- 
ter sudden removal of an electron from the HOMO of 
some organic unsaturated nitroso compounds known to 
possess low-lying satellites. Our results show that in all 
cases the initially created hole charge remains stationary 
but the system reacts by an ultrafast cyclic excitation- 
deexcitation process which leads to severe changes of 
the charge. In the presented example, the molecule 2- 
Nitroso[l,3]oxazolo[5,4-d][l,3]oxazole, the ir — ir* excita- 
tion following the removal of the HOMO electron takes 
place on a sub-femtosecond time scale, the period of the 
excitation-deexcitation alternations being about 1.4 fs. 
In real space the processes of excitation and de-excitation 
represent ultrafast derealization and localization of the 
charge. In other systems with more localized 7r* or a* 
orbitals one could anticipate by the same mechanism di- 
rected ultrafast charge migration from one site to an- 
other. 

At the end we would like to note that this dynamical 
behavior of the electronic cloud can, of course, influence 
the nuclear dynamics which will come into play at later 
times since, at least within the Born-Oppenheimer ap- 
proximation, the electronic motion governs the effective 
potential seen by the nuclei. 
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